CLNS 11/2079 
CLEO 11-05 

Branching fractions for T(3S') — )- TT^ht and ?/?(2S') — )- tt^/Ic 

J. Y. Ge/ D. H. Miller/ I. P. J. Shipsey/ B. Xin/ G. S. Adams,^ J. Napolitano,^ 
K. M. Ecklund,^ J. Insler,^ H. Muramatsu,^ C. S. Park,'* L. J. Pearson,* E. H. Thorndike,* 
S. Ricciardi,^ C. Thomas,^' ^ M. Artuso,^ S. Blusk,^ R. Mountain,'' T. Skwarnicki,^ 
S. Stone,'^ L. M. Zhang,^ G. Bonvicini,^ D. Cinabro,*^ A. Lincoln,^ M. J. Smith,^ 
P. Zhou,s J. Zliu,s P. Naik,9 J. Rademacker,^ D. M. Asner,*°'* K. W. Edwards,i° 
K. Randrianarivony,^° G. Tatishvili,*°' * R. A. Briere,^* H. Vogel,*^ P. U. E. Onyisi,^^ 
J. L. Rosner,*2 J. p. Alexander,^^ q q Cassel,*^ S. Das,^^ Ehrlich,*^ L. Gibbons,*^ 
S. W. Gray,*3 D. L. Hartill,*^^ B. K. Heltsley,*^ D. L. Kreinick,!^ V. E. Kuznetsov,*^ 
J. R. Patterson,*^ D. Peterson,*^ D. Riley,*^ A. Ryd,^^ j Sadoff,!^ x. Shi,*^ 
W. M. Sun,^3 j_ Yelton,!* P. Rubin,^^ N. Lowrey,!^ S. Melirabyan,!^ M. Selen,^^ 
J. Wiss,i6 J. Libby,*^ M. Kornicer,*^ R. E. Mitchells C. M. Tarbert,*^ D. Besson,*^ 

T. K. Pedlar,20 D. Cronin-Hennessy,^! J. Hietala,^! S. Dobbs,^^ Z. Metreveli,^^ 
K. K. Seth,22 A. Tomaradze,^^ T. Xiao,^^ L. Martin,^ A. Powell,^ and G. Wilkinson^ 

(CLEO Collaboration) 

^Purdue University, West Lafayette, Indiana J^lQOl, USA 
'^Rensselaer Polytechnic Institute, Troy, New York 12180, USA 
^Rice University, Houston, Texas 77005, USA 
"^University of Rochester, Rochester, New York 14627, USA 

^STFC Rutherford Appleton Laboratory, 
Chilton, Didcot, Oxfordshire, 0X11 OQX, UK 
^University of Oxford, Oxford 0X1 3RH, UK 
"^Syracuse University, Syracuse, New York 13244, USA 
^ Wayne State University, Detroit, Michigan 48202, USA 
^University of Bristol, Bristol BS8 ITL, UK 
Carleton University, Ottawa, Ontario, Canada KIS 5B6 
Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA 
^"^ University of Chicago, Chicago, Illinois 60637, USA 
Cornell University, Ithaca, New York 14853, USA 
University of Florida, Cainesville, Florida 32611, USA 
^^Ceorge Mason University, Fairfax, Virginia 22030, USA 
University of Illinois, Urbana- Champaign, Illinois 61801, USA 

^'^ Indian Institute of Technology Madras, 
Chennai, Tamil Nadu 600036, India 

^^Indiana University, Bloomington, Indiana 47405, USA 

University of Kansas, Lawrence, Kansas 66045, USA 

'^^ Luther College, Decorah, Iowa 52101, USA 

"^^ University of Minnesota, Minneapolis, Minnesota 55455, USA 

'^'^Northwestern University, Evanston, Illinois 60208, USA 

(Dated: June 17, 2011) 



1 



Abstract 

Using e+e^ collision data corresponding to 5.88 x 10^ T(3S') [25.9 x 10^ ^(25*)] decays and 
acquired by the CLEO III [CLEO-c] detectors operating at the Cornell Electron Storage Ring, 
we study the single-pion transitions from T{3S) [ip(2S)] to the respective spin-singlet states hh^c]- 
Utilizing only the momentum of suitably selected transition-7r'^ candidates, we obtain the upper 
limit i3(T(35) TT^h) < 1.2 x lO'^ at 90% confidence level, and measure B{il){2S) n'^hc) = 
(9.0 ± 1.5 ± 1.3) X 10~^. Signal sensitivities are enhanced by excluding very asymmetric 7r° — > 77 
candidates. 



* Now at: Pacific Northwest National Laboratory, Richland, WA 99352 
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Hadronic transitions between heavy quarkonium states provide a rich tableau of oppor- 
tunities for experimental and theoretical investigations [1]. Such transitions can supply im- 
portant production mechanisms for the lower-lying state, thereby allowing measurement of 
its mass, width, and decay channels, and can also enable nonperturbative QCD calculations 
of transition rates to confront experiment. In particular, observation and exploration of the 
bb [cc] spin-singlet states hb[c] have depended strongly upon hadronic transitions from vector 
quarkonia produced in e+e" collisions. The /if,[c] states^ were both expected and measured 
to have masses near the spin-weighted averages of the respective spin-triplet states Xb[c]j- 
Observations of 'ip{2S) — ?■ n^hc were reported first by CLEO [2] and later by BESIII [3]. 
Evidence for T{3S) — ?■ n^hf, has been reported by BABAR [4]. Charged dipion transitions 
from higher-lying vector states were first shown to produce he by CLEO [5] and, later, both 
/lb (IP) and h{2P) by Belle [6]. 

In this article we report on attempts to measure the branching fractions for the isospin- 
violating single- 7r° transitions from T(3S) [ipi^S)] to hb[c]- Previously, BABAR [4] measured 
the product branching fraction B{T{3S) -> n'^h) x B{hb 7r/b) = (4.3 ± 1.1 ± 0.9) x 10"^ 
finding evidence for a signal with significance of 3.3a, consistent with predictions B(T{3S) 
n^h) ^ (1.0-1.6) X 10-3 [7, 8] and B{hb ^rVb) = 41% [9]. Meanwhile, BESIII [3] measured 
B{ijj{2S) — )■ Ti^hc) = (8.4±1.3±1.0) x 10"^ (in conjunction with a simultaneously determined 
B{hc — )■ jrjc) = (54.3±6.7±5.2)%), compatible with both potential model estimates [10] and 
a non-relativistic effective field theory prediction [11] (see also [12]) which found charmed 
meson loop contributions to mass shifts to be atypically small. 

The two data sets used in this work were collected in e~^e~ collisions at the Cornell 
Electron Storage Ring, at the center-of-mass energies of the T(3S') and ip{2S) resonances, 
corresponding to (5.88±0.12) x 10^ T{3S) decays [13] and (25.9±0.5) x 10^ ^Ij{2S) decays [14]. 
Events were recorded using the CLEO III and CLEO-c detectors for T(35') and 'ipi^S) 
datasets, respectively. Both configurations are equipped with an electromagnetic calorimeter 
consisting of 7784 thallium-doped cesium iodide crystals and covering 93% of solid angle, 
initially installed in the CLEO II [15] detector configuration. The energy resolution of 
the crystal calorimeter is 5% (2.2%) for 0.1 GeV (1 GeV) photons. Calorimeter angular 
resolution is ~10 (5) mrad at = 100 MeV (1 GeV), and does not significantly contribute 
to 77 mass resolution for the soft vr*^ candidates considered in this analysis. The CLEO III 
tracking system [16] consists of a silicon strip vertex detector and a large drift chamber; a 
six-layer wire vertex detector replaced the silicon in the CLEO-c configuration [17]. The 
trackers achieve charged particle momentum resolutions of 0.35% and 0.6% at 1 GeV/c in 
1.5 T and 1.0 T axial magnetic fields, respectively. 

In both studies we demand the presence of only the transition pion (via — )■ 77) in 
hadronic events, and impose no restriction on the hb[c] decay other than the global event 
selections described later. Branching fractions for hb[c] — )■ 'JVbic] are known to be large, so 
that hb[c] searches can and do reasonably seek to require either the radiative photon and/or 
a reconstructed rib[c]] we avoid such an approach here. Instead, the magnitude of the hb[c] 
signal is inferred from the size of any enhancement in the distribution of Mrec(vr°), the mass 
recoiling against the putative transition 7r° (Mj.ec(7r°) = a/ {pres — Ptto)'^ where pres and p^o 
are the initial T{3S)[ip{2S)] and the vr" four-momenta, respectively). The main challenge is 
to design selection criteria to simultaneously preserve signal strength while suppressing the 
unavoidably large backgrounds. 



Throughout this article, lowest radial excitations are implied unless explicitly indicated. 
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TABLE I. Four key values used in each MC-based optimization of selection criteria and some 
related energies of interest. 



Item 


Units 


T(3^) 


^(2^) 




MeV/c2 


9900.0 


3525.28 




MeV/c2 





0.86 


f3{h[c] 7^fe[c]) 


% 


38.0 


37.7 




10-4 


16 


8.4 




MeV 


446 


159 




MeV 


433, 452 


128, 171 



We determine our selection criteria and fit procedure based on Monte Carlo (MC) simu- 
lations of resonance decays, continuum (e^e" — )■ 7* — )■ qq) background, and an off-resonance 
sample, 20.7 pb^^ taken ~16 MeV/c^ below the ^(25*), so as to avoid bias. The MC gen- 
eration utilizes EvtGen [18], the values shown in Table I, and the most recent relevant 
branching fraction, mass, and width measurements [13, 19, 20] and predictions (where mea- 
surements do not exist), followed by a GEANT-based simulation [21] of the two detector 
configurations. First, global event selections for T{3S) and ipi^S) are performed (with 92% 
and 95% efficiencies) as in Refs. [13] and [22], respectively. In order to isolate reliably- 
triggered resonance decays and suppress e~^e~ , , and r+r^ final states, the selections 
demand that there be well-reconstructed charged particle tracks and that the total mea- 
sured energy in each event be consistent with that of the initial e~^e~ state. The copious 
background from ip{2S) — )■ -k^-k^ J /tp is reduced by requiring that events have no oppositely- 
charged pair of particles with dipion recoil mass near that of J/ip. 

We reconstruct the transition 7r° — i- 77 candidates based on pairs of showers in the 
calorimeter that are not matched to the projected trajectory of any charged particle. Showers 
must be located well within the boundaries of the crystal calorimeter barrel (| cos^| < 0.81) 
for 'iIj{2S) -> Ti^hc, or barrel and endcaps (0.85 < | cos 6^1 < 0.93) for T(3S') — > vr^/ift, 
where 9 is the polar angle with respect to the positron beam direction. Also, the showers 
must have > 30(50) MeV in the barrel (endcap) calorimeter. To reduce background 
from non-photon hadron-induced showers, photon candidates are also required to have a 
lateral shower profile consistent with that of an isolated electromagnetic shower. Defining 
the 7r° mass "pull" (= {M^^ — M^o)/AM^^, where AM^^ is the photon-energy-dependent 
resolution on 77 invariant mass, typically ~5-7 MeV/c^), we find that restricting the mass 
pull to [—3.0, +2.5] optimizes sensitivity to signal. The asymmetric mass pull range accounts 
for the presence of a low-side tail in M^^ caused by lateral and longitudinal shower leakage 
from the crystals assigned to the photon candidate's shower. If a daughter photon is shared 
with more than one vr" candidate, the pair with smaller mass pull is chosen. We then 
kinematically constrain M^^ to the known M^o to improve 7r° momentum resolution. 

By the very nature of such an inclusive measurement, most of our selected events will 
be background, and any inference of an hh^c] signal depends strongly on two characteristics 
of the background: first, that it has smooth Mrec(vr°) dependence in the vicinity of the hb[c] 
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FIG. 1. Distributions of the -k^ decay angle a for MC (or data, where noted below) samples 
scaled to the size of our datasets: (a) T(3S') — vr'^/if, signal (dashed), all T(35) MC decays except 
signal (solid), T(35) -yxbj (dotted), and e+e" ^ 7* gg (dash-dot), for 9875 < Mrec(vr°) < 
9925 MeV/c^; (b) ^(25) ^ 7r°/ic signal (dashed), ah ^(25") decays except signal (solid), V'(2S') 
7Xcj (dotted), and below-V'(2S') continuum data (dash-dot), for 3520 < Mrec(7r°) < 3530 MeV/c^. 
Vertical lines and arrows show selected regions of | cosa|. 



mass, so that extrapolation of its shape underneath any signal can be made with confidence, 
and second, that its magnitude can be reduced enough to observe a peak of adequate statis- 
tical significance. Enormous combinatoric smooth backgrounds are present from tt*^ — >■ 77 
and rj — 77, which can arise from e^e~ — )■ qq or at various stages of resonance transitions or 
decays. Although without structure in Mrec(7r°), these must be suppressed without sacrific- 
ing too much signal. Moreover, significant non-smooth structures in Mrec(7r°) are necessarily 
present as well. Because the spin-singlet masses are near the spin-weighted averages of the 
respective spin-triplet states, the transition-vr*^ energy is always close to the photon energies 
from the electric-dipole (El) transitions T{3S)['iIj{2S)] 1Xb[c]j, as shown in Table 1. An 
El photon will frequently be paired with one of the multitude of very low energy photons 
in most events to form a mass close to that of a true 7r°. As the above-mentioned El pho- 
tons are monochromatic, non-smooth contributions in Mrec(vr°) can arise near the h^c] mass. 
These sources of fake 7r° candidates make extrapolating a reliably known background shape 
underneath any hbic] signal systematically challenging unless mitigating measures are taken. 

Our chosen method for fake-TT*^ background suppression is to restrict the values of the 
7r° — > 77 decay angle a, taken as the angle in the tt" center-of-mass frame between either 
photon and the vr^-boost direction. True 7r° decays have a uniform distribution in | cosa|. 
Values of |cosq;| near unity imply an asymmetric decay, with one of the photons being 
very soft; it is these candidates which give rise to most of the background. MC studies 
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FIG. 2. Solid histograms show distributions in Mj.cc{tt^) for MC samples of signal tt^'/i^j^j transitions 
at modeled levels plus (a) jXbj, and (b) 7Xc2 backgrounds, scaled to the size of our datasets, with 
our restrictions on decay angle [(a) |cosa| < 0.7, (b) < 0.5] (lower pairs of histograms) and 
with I cosq| < 1.0 (upper pairs). Dashed histograms show the contribution of (a) 'jXbj, and (b) 
7Xc2-only events for both selections on | cosq|. 



indicate that | cosa| < 0.7 [0.5] provides the best compromise between sensitivity and back- 
ground rejection for T{3S) — > ir^hb [4'{2S) — )■ vr^/ij. The tighter value for ip(2S) — )■ n^hc 
reflects the order-of-magnitude larger values of B{ip{2S) '-yXcj) l^^ipi^S) — )■ 7r°/ic) rela- 
tive to i3(T(3S') — )■ 7Xfej)/'B(T(3S') — > 7r°/i{,) and the consequent need to suppress the 7Xcj 
backgrounds more severely. All these points are demonstrated in Figs. 1 and 2, where it 
can be seen that backgrounds congregate at larger | cosq;| and cause structure in Mrec(7r°) 
unless suppressed. The signal histograms in Fig. 1 fall with increasing | cosa| due to the 
requirement that > 30(50) MeV in the barrel (endcap) calorimeter, the Mj-edir'^) range re- 
striction for these plots, and other event selection criteria. Exclusion of asymmetric 7r° — )■ 77 
decays also has the advantage of improving Mrec(vr°) resolution (which enhances sensitivity) 
because such decays include softer photons, which have poor relative energy resolution. 

Signal extraction is accomplished by fltting each Mi.ec(7r°) distribution with suitable bin- 
ning and range for the combination of smooth background and peaking signal component 
shapes with flxed M{hi,[c]) and floating normalizations for each; choices of mass, binning, 
range, and shapes appear in Table II. The hb signal shape is found to be best represented by 
a reversed Crystal Ball line (CBL) shape (an ordinary CBL shape [23] with the power-law 
tail on the high side) while flxing the shape parameters based on signal-only flts to MC 
samples. A third-order polynomial for the hb background flts the data well. The he shape 
is chosen instead as a double Gaussian with independent means (the mean of the broader 
Gaussian is shifted higher) based on studies of our signal MC sample. The different signal 
shapes reflect those expected for very slow {he) and faster {hb) 7r° mesons, and reflect that 
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TABLE II. Features of the fits for signal extraction. "Order" in background shapes refers to order 
of polynomials. See text for other details. 



T(3^) i){2S) 



Item 

^reci^^^) binning 
^rec(7r^) fit range 
Background shape 
M(/i,[,]) for fit 
Signal shape 



4 MeV/c^ 
9.8-10.0 GeV/c^ 
S"^*^- order 
9900 MeV/c2 
Reversed CBL 



1 MeV/c^ 
3505-3545 MeV/c^ 
ARGUS+2^'i-order 
3525.42 MeV/c^ 
Double Gaussian 



the calorimeter resolution is more symmetric at lower energies and develops a low-side tail 
at higher energies due to shower leakage; the tail moves to the high side in recoil mass. Since 
Et^o from ip{2S) — )■ vr^/ic is close to the kinematic limit of M{7t^), the he background shape 
is represented by an ARGUS function [24] , which effectively models reduced phase space as 
a function of increasing mass, plus a second-order polynomial. Our MC studies show we can 
extract branching fractions consistent with what we input to the MC samples with these fit 
procedures. 

Results from the fits appear in Figs. 3 and 4 and Table 111. No signal is seen for 
T(35') — )■ n^hb, and the quoted upper limit integrates over physical branching fractions 
only and includes systematic effects (see below). An unambiguous signal is observed for 
ipi^S) — TT^hc- The data points above the fit level near 3513 MeV/c^ in Fig. 4(b) have a 
width narrower than the detector resolution, and hence must constitute a statistical fluctua- 
tion. The statistical significances shown are computed as In {L„o/L„), where L^o and 
are likelihood values from fits of Mrec(7r'') without and with signal shape components, 
respectively. If the he mass is allowed to float, 3525.9 ±0.3 MeV/c^(a mass value l.la larger 
than the world average [19]) is obtained which results in negligible change in fitted yields 
with respect to the case when the mass is fixed to the world average. The final branching 
fractions are obtained as = Nevt/i^N^es) where Nevt is the number of signal events ex- 
tracted from the fit, iVres is the number [13, 14] of resonance decays in the dataset, and e is 
the reconstruction efficiency obtained from performing similar fits on MC samples. 

We consider a variety of sources for systematic uncertainties on the branching fractions 
obtained from the fits, and summarize them in Table IV along with estimates of their con- 
tributions. The general approach is to vary the important selection criteria or fitting choices 
over reasonable ranges and note any resulting variations in B{T{3S)[ip{2S)] vr^/ib^c]) be- 
yond expected statistical changes, and then to add all such effects in quadrature. The 
dominant systematic effects are quite different for h^ and he- For hb, the source that stands 
out is the Mrec(vr°) fit range. For h^ the largest contributions come from the functional form 
of the background shape, uncertainty in T{hc), and our understanding of vr*^ resolution in 
data and MC simulation. As Mj-edn^) resolutions are larger than the hb[c] mass uncertain- 
ties [3, 4, 6, 19], there is no need for separate systematic errors from such variation. 

For the h Mrec(7r°) fit range, we consider alternate ranges 50 MeV/c^ wider and nar- 
rower, symmetrically around our chosen M{hb), observing excursions as large as noted in 
Table IV. For background shape, a fourth-order polynomial is tried instead of a third-order; 
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FIG. 3. (a) Fit to Mrec(^°) for T(3S) ^ TT^hj, for fixed M{hb) = 9900.0 MeV/c^. The value 
from this fit is 26.4 for 50 data points (minus 5 parameters) with confidence level of 98.8%. (b) 
The fitted background-subtracted spectrum (solid curve). The dashed curve corresponds to the 
upper limit on signal candidates at 90% CL (< 1439 events or ^(T(35) 7r°/ib) < 11 x 10""^ at 
CL). 



for signal shape, we allow a double Gaussian instead of the reversed CBL shape (effects 
of imperfectly understood Mi.ec(7i"'^) resolution are also accounted for in this variation). To 
test the dependency of our result upon the predicted resolution in Mi.ec(7r°), we decrease the 
smearing predicted by the MC by 8.5% less than predicted (which no longer gives reasonable 
agreement between data and MC samples), and consider half of the change in its measured 
branching fraction as a possible systematical bias. For binning of Mrec(7r''), we vary from 
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FIG. 4. (a) Fit to the ip{2S) — t- ir^hc Mrec(7r*^) data with (sohd curve) and without (dashed) a 
signal. The value from this fit is 41.0 for 40 data points (minus 5 parameters) with confidence 
level of 22.5%. (b) As in (a) but with the background fit from (a) subtracted. 

4 to both 2 and 6 MeV/c^ for hf,. We allow B{hb — )■ •yrib) to vary from 0% and up to 
100% because its size is unknown and has a small but nonzero effect on photon multiplicity 
and therefore upon the efficiency of the signal vr" reconstruction. We also account for the 
uncertainty in iVj-es- 

For the ip{2S) — )■ n^hc fit, the background shape is alternately fit to either first- or 
third-order polynomials instead of the second-order. The he signal size appears to have an 
approximately linear dependence on the assumed T{hc), behaving as B{ip{2S) — >■ n^hc) = 
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TABLE III. Signal-extraction efficiencies and final measured event yields and branching fractions, 
the latter including systematic uncertainties. The upper limit integrates over physical values only. 
See text for details. 



Item T{3S) 7/>(25) 

e(7r*^/i5[c]) 21.3% 12.6% 

TVevt 139 ±821 2943 ±501 

Significance 0.2cr 5.9cr 

A^res(lO^) 5.88 ±0.12 25.9 ±0.5 



^(7r°/i6[c]) (10-^) < 12 at 90% CL 9.0 ± 1.5 ± 1.3 



TABLE IV. Summary of relative systematic uncertainties in percent on B{T{3S) — )■ ir^hi,) and 
B{^{2S) — )• TT^hc). Entries marked "• • •" make negligible contributions. 



Source 


hb 


he 


Background shape 


0.8 


9.3 






7.8 


Fit range 


19.7 




Binning 


10.9 




Signal shape 


1.7 




TT^ resolution 


2.0 


6.6 


A^res 


2.0 


2.0 


^{h[c\ lVb[c\) 


2.5 


4.1 


Efficiency (MC statistics) 


0.4 


0.6 



I cosal 

Quadrature sum 22.9 14.6 



[7.6 + lAT{hc)/To] X 10"'^, where Tq = 0.86 MeV/c^ is the chosen width. We then vary the 
width by ±50% of 0.86 MeV/c^ to estimate a systematic error. We account for uncertainty 
in calorimeter resolution by varying it over ranges that still represent the data reasonably 
well, as in our study. We vary | cosalmax (0.5 ± 0.1), the fit range (3505-3550 MeV/c^; 
we could not go lower than 3505 MeV/c^ due to the contamination from '0(25') — )■ 7Xc2), 
suppression of events from 'ip{2S) — t- tt^ti^ J/tp, bin width (1.0 ± 0.5 MeV/c^), as well as 
doubling and halving the assumed B{hc — )■ 7^7c), which includes the value measured by 
BESIII [3]. 

In conclusion, we have measured branching fractions for T{3S)[iIj{2S)] — ir^h^c] as shown 
in Table III. The hh upper limit is dominated by statistical uncertainties and supersedes the 
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previous CLEO limit, B[T{3S) n^h] < 2.7 x lO'^ at 90% CL [25]. If we combine 
the product branching fraction from BABAR [4] with our resuh, considering only physical 
values of branching fractions, we infer B{hb — )■ 7?7b) > 24% at 90% CL, consistent with 
predictions [9]. The he result is consistent with the value from BESIII [3]. 
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